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Abstract 

Using e+e" -> D~D*+ and D*~D+ interactions at 4170 MeV collected with the CLEO- 
c detector, we investigate the semileptonic decays Df — > /o(980)e + zv, and — > ^e + i^. 
By examining the decay rates as functions of the four-momentum transfer squared, q 2 , 
we measure the ratio \$(D+ -► / (980)e + i/)£(/ tt+tt-)] / [|l(£>+ -► <f>e+v)B(<l> -► K+K~j 
at g 2 of zero to be (42±11)%. This ratio has been predicted to equal the rate ratio 
[B(B S ^ J/ijfo)B(f ^ir+ir-)}/[B(B s ^ J/i>(P)B(<j)^ K+K-)], thus indicating that the CP 
eigenstate J/ipfo could be useful for measuring CP violation via B s mixing. Assuming a simple pole 
model for the form factor |/+(<? 2 )| in the foe + v decay, we find a pole mass of (1.7l '|±0.2) GeV. We 
also determine the /o mass and width as (977±J 1 ±1) MeV, and (9lj^±3) MeV, respectively. In ad- 
dition, we present updated results for B{D+ -► / (980)e + z/)£(/ -► vr+vr-) = (0.20±0.03±0.01)%, 
and B(D+ — ► c/>e + ^) = (2.36 ± 0.23 ± 0.13)%. Assuming that the /o wavefunction is a com- 
bination of strange and non-strange quark-antiquark components, we use our measurement for 
B(Df — ► /o(980)e + z/) to extract a value of the mixing angle that we find consistent with | ss) 
dominance, adding to the mystery as to why the /o decays predominantly to two pions rather than 
two kaons. 

PACS numbers: 13.20.Fc, 12.38.Qk, 14.40.Lb 
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I. INTRODUCTION 



In this article we present a study of semileptonic decay of the meson into /o(980)e + z/ 
and also to (f>e + v. The /o(980) meson, a scalar, though well established experimentally, 
has a relatively uncertain mass and width [1]; in addition there have been claims that the 
quark content may be a mixture of a traditional quark- ant iquark with a four-quark system 
[2]. Semileptonic decays provide a pristine environment where the fo is produced by an 
isoscalar combination of s and s quarks. Evidence for / (980) in semileptonic decays was seen 
by BaBar via the /o(980) — > K + K~ channel, where interference was observed by an S-wave 
with the dominant P-wave decay [3]. The first measurement of the Df — > /o(980)e + z/, 
fo — > 7t + tt~ product branching fraction was recently made by CLEO [4]. The semileptonic 
decay diagram is shown in Fig. 1. Here we investigate this mode using a data sample of 
600 pb _1 , approximately double the original size. This larger sample is sufficient to allow 
us to determine the fo mass and width and measure the semileptonic decay form factor as a 
function of the invariant four- momentum transfer squared, q 2 , between the D s and the fo- 
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FIG. 1: The Feynman diagram for semileptonic Df decay into a or /o(980) meson. 

CP violation measurements in the B s system have concentrated on the final state J/i/xf), 
with — > K + K~ . For a review see Ref. [5]. Since this mode is not a CP eigenstate an 
angular analysis is required to separate the CP even and CP odd parts, and measure the CP 
violating phase —2/3 s . 

In what follows we will use a simple notation for the decay width of a particle multiplied 
by the branching fraction for the decay of one of its daughters. For example, 

T(D+ - 7T + /0, fo -> 7T + 7T-) EE T(D+ - 7r + /o(980))£(/ (980) -> 7T + 7T-). (1) 

Stone and Zhang [6] have suggested that the J/ip4> mode may also contain an S-wave 
K + K~ system at the mass with a rate that could be as large as ~5% that of the 0. This 
S-wave, if it should be significant, would require additional parameters in the fit to extract 
— 2(3 S . They also suggest that the final state J/ijifo may be a useful alternative; since it is 
a CP-eigenstate angular analysis is not required. Stone and Zhang estimated the branching 
ratio B s — > J/ipfo assuming equality of the ratios 

T(B° s ^J/^fo, /q^tt+tt-) _ T(D+^f n+, f -> tt+tt-) ^ 

K m - r ( £o ^ W) ^ K+K-) ~ r {D + ^ 07r + ^ K+K -) - ^ 6[) > /o ■ w 

The phase spaces for the _B° and decays are quite similar, but the spin structure is not. 
In these Df decays both the and the fo are produced opposite a spin-0 pion, rather than 
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a spin-1 J/ip. This consideration prompted Stone and Zhang to suggest that measuring 
the ratio of decay widths of semileptonic -D+ decays containing either an f or a <fi at four- 
momentum transfers, q 2 , equal to zero would give a superior prediction [6]. Specifically they 
propose that 

^(Df -> / (980)e+z/, f -> 7T+7T-) | 92=0 

The point q 2 equal to zero is chosen to maximize the allowed phase space in order to make 
it as close as possible to that available in B° s — > J/ip(fi (or / ) decay. In this paper we will 
present measurements of this ratio, the form factor in the — > /o e+z/ channel, the / mass 
and width, and update the previously published CLEO branching fractions for these two 
semileptonic decay modes [4]. 



II. EXPERIMENTAL METHOD 
A. Selection of D s Candidates 

The CLEO-c detector [7] is equipped to measure the momenta and directions of charged 
particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH) [8], 
detect photons and determine their directions and energies. 

In this study we use 600 pb™ 1 of data produced in e + e~ collisions using the Cornell Elec- 
tron Storage Ring (CESR) and recorded near a center-of-mass energy (Eqm) of 4.170 GeV. 
At this energy the e + e~ annihilation cross-section into D~D* S + plus D* S ~D+ is approximately 
1 nb [9]. 

In this analysis we fully reconstruct a sample of D~ in several "tag" modes and then find 
candidate semileptonic decays in this sample. Mention of any specific decay implies the use 
of its charge-conjugate as well. The tag selection is identical to that used in our — > fi + u 
paper, that can be consulted for details [10]. Briefly, we select candidates on the basis of 
their beam-constrained invariant mass. Then we detect an additional photon candidate from 
the D* decay, and construct the missing mass squared, MM* 2 recoiling against the photon 
and the D~ tag 

MM* 2 = (E CM - E Ds - E,f - (p CM - p Ds ~ P 7 ) 2 , (4) 

where -Ecm (Pcm) is the center-of-mass energy (momentum), Ejj s (Pd s ) is the energy (mo- 
mentum) of the fully reconstructed D~ tag, and E 1 (p 7 ) is the energy (momentum) of the 
additional photon. In performing this calculation we use a kinematic fit that constrains the 
decay products of the D~ to the known D s mass and conserves overall momentum and en- 
ergy. All photon candidates in the event are used, except for those that are decay products 
of the D~ tag candidate. Regardless of whether or not the photon forms a D* with the tag, 
for real D*D S events MM* 2 should peak at the mass-squared. 

We list the number of signal events in each mode in Table I by finding the number of 
events within ±17.5 MeV of the D s mass. For ease of further analysis we sum all tag modes 
together, as shown in Fig. 2(a). 

The MM* 2 distributions for events in the D~ invariant mass signal region (±17.5 MeV 
from the D s mass) are shown in Fig. 2(b). In order to find the number of tags used for 
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TABLE I: Tagging modes and numbers of signal and background events, within ±17.5 MeV of the 
D~ mass for each mode, determined from two-Gaussian function fits to the invariant mass plots, 
and the number of tags in each mode including the 7 from the D* — > jD s transition, within an 
interval 3.782 < MM* 2 < 4.0 GeV 2 , as determined from fits of the MM* 2 distributions (see text) 
to a signal Crystal Ball function (see text) and two 5th order Chebychev background polynomial 
functions. 



Mode 


Invariant Mass MM* 2 




Signal Background Signal Background 



K+K-JT- 


26534±274 


25122 16087±373 


39563 


K S K- 


6383±121 


3501 


4215±228 


6297 


i]tt~ ; i] — ► 77 


2993 ± 156 


5050 


2005±145 


5016 


ij'tt - ; rj — > ir + ir~7], 77 — > 77 


2293±82 


531 


1647±131 


1565 




11649±754 


78588 


6441±471 


89284 


7T + 7T~7r~ 


7374±303 


60321 


5014±402 


43286 


K *- K *o. K *- Kgw~, K*° -»• K + tt- 


4037±160 


10568 


2352±176 


12088 


vp~;v^ 77, P~ n-TT 


5700±281 


24444 


3295±425 


24114 


rfn~; if -> p°7, 


3551±202 


19841 


2802±227 


17006 


Sum 


70514 ± 963 


227966 43859±936 


238218 



further analysis we perform a two-dimensional binned maximum likelihood fit of the MM* 2 
distribution and the invariant mass distribution in the interval ±60 MeV from the D s mass 
and 3.5 < MM* 2 < 4.25 GeV 2 . The background has two components, both described by 
5th order Chebyshev polynomials in MM* 2 ; the first comes from the background under the 
invariant mass peak, defined by the sidebands, and the second is due to multiple photon 
combinations. In both cases we allow the parameters to float. 

We find a total of 43859±936±877 events within the interval 3.782 < MM* 2 < 4.0 GeV 2 
and having an invariant mass within ±17.5 MeV of the D s mass, where the total number 
of events is the sum of the yields from the fits to each mode as shown in Table I. The first 
uncertainty in the total is statistical and the second is systematic. 

B. Signal Reconstruction 

We next describe the reconstruction of Df — > /o e+z/ ) fo — * vr + 7r~, and also — > 
0e + z/, — > K + K~ . We select events within the MM* 2 region shown in Fig. 2(b) for further 
analysis. We note that the limits are rather wide. We use this selection because the back- 
ground in the signal side is rather small and the errors are minimized by taking as many 
tags as possible. 

Candidate events are selected that contain a charged track of opposite sign to the D~ tag 
that is positively identified as a positron. Electrons and positrons are identified on the basis 
of a likelihood ratio constructed from three inputs: the ratio between the energy deposited in 
the calorimeter and the momentum measured in the tracking system, the specific ionization 
dE/dx measured in the drift chamber, and RICH information [11]. Our selection efficiency 
averages 0.95 in the momentum region 0.3-1.0 GeV, and 0.71 in the region 0.2-0.3 GeV. The 
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FIG. 2: (a) Invariant mass of D~ candidates summed over all decay modes and fit to a two-Gaussian 
signal shape plus a straight line for the background. The vertical dot-dashed lines indicate the ±17.5 
MeV definition of the signal region, (b) The MM* 2 distribution summed over all modes. The curves 
are fits to the number of signal events using the Crystal Ball function and two 5th order Chebychev 
background functions; the dashed curve shows the background from fake D~ tags, while the dotted 
curve in (b) shows the sum of the backgrounds from multiple photon combinations and fake D~ 
tags. The vertical dashed lines show the region of events selected for further analysis. 



average fractions of charged pions and kaons incorrectly identified as positrons averaged over 
the relevant momentum range are approximately 0.1%. We also require an additional pair 
of tracks with opposite charge that are both identified as pions or kaons using the dE/dx 
and RICH systems. 

Since we are searching for events containing a single missing neutrino the missing mass 
squared for the f mode, MM 2 , evaluated by taking into account the observed e + , 7r + , n~ , 
D~ , and 7 should peak at zero; the MM 2 is computed as 



MM 2 = (E C m - E Ds — E 1 — E e — E n + - E n -f - (p C M - Pd s - P 7 - Pe - Pn+ ~ Pn-f , 

(5) 

where E e (p e ) are the energy (momentum) of the candidate positron, E n (p n ) are the energy 
(momenta) of the candidate pions, and all other variables are the same as defined in Eq. (4). 
A similar equation applies for the mode with the pions replaced by kaons. 

We also make use of a set of kinematical constraints and fit each event to two hypotheses, 
one of which is that the D~ tag is the daughter of a D*~ , and the other that the D* + 
decays into with the subsequently decaying into either 7r + 7r~e + z/, or K + K~e + v. 
The kinematical constraints, in the e + e~ center-of-mass frame, are 

Pd s + Pd* s = 0, (6) 
Eqm = E Ds + E D * , 

g» Ml, - Ml. m gcM _ - Ml, 

M D * - M Ds = 143.8 MeV. 
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In addition, we constrain the invariant mass of the D~ tag to the known D s mass. This gives 
us a total of seven constraints. The missing neutrino four-vector needs to be determined, so 
we are left with a three-constraint fit. We perform an iterative fit that minimizes x 2 ■ As 
we do not want to be subject to systematic uncertainties that depend on understanding the 
absolute scale of the errors, we do not make a x 2 cut but simply choose the photon and the 
decay sequence in each event with the minimum x 2 ■ 

We model the MM 2 signal distributions for both /o e+z/ an d 0e + z/ final states as the sum 
of two Crystal Ball (CB) functions centered at zero [12] plus a wide Gaussian shape that 
serves to model the tails. The Monte Carlo simulations are shown in Fig. 3. The results are 
summarized in Table II. 
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FIG. 3: MM 2 distributions from Monte Carlo simulation for (a) Df — > foe + v and (b) Df - 
The curves are the sum of two CB functions with means fixed at zero, n fixed at 1.8 and a values 
constrained to be equal, and a single Gaussian shape. (The parameters are listed in Table II.) 



TABLE II: Parameters of signal fits to two CB functions plus a Gaussian shape; "J 7 " indicates the 
fraction of the yield in each fit component. The parameter n used in both CB functions is fixed at 
1.8. The parameter a is kept the same for both CB functions. The mean and width of the Gaussian 
is kept the same for both decay modes. All means and cr's are in units of GeV 2 . 



1 st CB function 



a 



2 nd CB function 



Gaussian 
J r (%) mean a 



K+K- 



63 (6.30 ±0.15) x 10" 3 1.417±0.014 30 (1.67 ± 0.07) x 10" 2 
45 (5.61 ± 0.35) x 10~ 3 1.610±0.003 37 (1.32 ± 0.09) x 10~ 2 



7 -0.0253 0.0992 
18 -0.0253 0.0992 



We proceed by performing a simultaneous fit to the D~ invariant mass, using a mass 
range ±70 MeV from the nominal mass, and the MM 2 for 7r + 7r~e + z/ and K + K~ e + v to the 
Monte Carlo generated fitting functions letting only the normalizations float. 1 The resulting 



1 This is the same procedure as used in our D+ — ► fi + v analysis. 
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MM 2 distributions are shown in Fig. 4. Here and in subsequent analyses we require that the 
7r + 7r~ invariant mass is above 0.6 GeV, in order to eliminate the K^e + v channel, and limit 
the K + K~ invariant mass to be below 1.08 GeV. 

Our next step is to determine the range of MM 2 to select for further analysis. Thus, we 
consider the backgrounds in our MM 2 sample. In general the background arises from two 
sources: one from real decays and the other from the background under the single-tag 
signal peaks (fake D~). For 0e + z/ only the background from fake D~ is significant. For 
foe + u we show in Fig. 4(a) the background from the invariant mass sidebands, suitably 
scaled. We also find that there are two sources of real -D+ background, namely r/e + z/ decays 
where the rf — > p7, and a very small component of 0e + z/, — > 7r + 7r~7r° when the tt° is 
ignored. We also show these backgrounds as a function of MM 2 . We use a second order 
polynomial to parameterize the fake D~ background, and the other two backgrounds are 
parameterized by the sum of two bifurcated Gaussian shapes determined by simulating the 
background processes. (A bifurcated Gaussian shape has a different widths below and above 
the mean.) We are able to fix the 7r + 7r~7r°e + z/ rate as 0.31%, using our previous measurement 
of -D+ — > <f)e + v, but allow the 7r + 7r~7e + z/ yield to float in the fit as our measured branching 
fraction is only accurate to ±35% [4]. 
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MM 2 (GeV 2 ) 

FIG. 4: (a) The MM 2 distribution for tt + tt e + v. The dotted curve shows the signal, the long- 
dashed distributions are from sidebands of the Dj candidate invariant mass distributions while the 
shorter dashed curve in (a) shows the rf — * 7r + 7r~7 background level, while the dashed-dot curve 
shows the (small) 7r + 7r _ 7r° background from (pe + v. The solid curve shows the total, (b) The MM 2 
distribution for K + K~e + v, where the total and the small side-band background is shown. 

For further analysis we restrict ourselves to the interval —0.04 < MM 2 < 0.04 GeV 2 , 
which is 88.4% efficient on foe + u and is 84.4% efficient on 0e + z/. In these regions we find 
42.9±6.7 7T + 7T~e + v signal events, 7.7±2.0 rf e + v events, 0.5 0e + z/ events (fixed), and 13.4±0.9 
fake D~ background. There are 107.0±9.8 K + K~e + v events of which 2.5±0.7 are fake D~ 
background. 

C. K + K~ and 7r+7r — Invariant Mass Distributions 

We next include the di-hadron mass as part of the joint fit to the D~ mass, and MM 2 . 
Figure 5 shows the 7r + 7r~ invariant mass distribution for 7i + 7i^e + u events, and the K + K~ 
invariant mass distribution in our selected sample of K + K~e + v events, within the interval 
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—0.04 < MM 2 < 0.04 GeV 2 . Both distributions are dominated by a single di-hadron reso- 
nance, either the /o(980) in the tt + tt~ mode or the <fi in the K + K~ mode. The backgrounds 
are included in the fit in the same manner as for the MM 2 distributions. We next find the 
size of the fo and signals. 
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1.08 



k + tz- Mass (GeV) K + K~ Mass (GeV) 

FIG. 5: Invariant mass distributions for (a) ir + ir~ and (b) K + K~ for semileptonic decays.. The 
dotted curve shows the signal, the long-dashed distributions are from sidebands of the D~ can- 
didate invariant mass distributions while the shorter dashed curve in (a) shows the rj' — ► 7r + 7r~7 
background level, and the dashed-dot curve in (a) shows the (small) 7r + 7r~7r background from 
(j)e + v. The solid curves show the totals. 



For the K + K~e + u, the signal fitting function consists of the sum of a P-wave Breit-Wigner 
shape whose mass and width are fixed to the values reported by the Particle Data Group 
(PDG) [1] convoluted with our Gaussian experimental resolution, and a fixed background 
function determined from the D~ sidebands described by the same <fr signal shape added to 
a wide component that is exponentially decreasing with mass. 

For the tt + tt~ mass distribution we include three separate shapes to fit the backgrounds. 
(1) The 7r + 7r~7 distribution is described by a bifurcated-Gaussian function; (2) the 7r + 7r~7r° 
distribution is described by a bifurcated-Gaussian shape plus a Breit-Wigner shape describing 
the p, here our knowledge of the (f>e + v decay rate allows us to fix the magnitude of this 
component; (3) the shape describing non-D~ background is an exponential plus a Breit- 
Wigner function describing the p, the size of this component is determined from the sidebands 
of the D~ signal. 

The mass and width of the /o(980) are not well determined. The PDG makes an estimate 
of the mass of (980±10) MeV, but does not explicitly average any measurements. They also 
state that the width is very model dependent and give a range of 40-100 MeV, again not 
forming an average [1]. Here we provide measurements in an extremely clean environment. 
We fit the 7r + 7r~ mass distribution, y/s = M wn with a Breit-Wigner type resonance form 
with the threshold behavior suggested in [13] 

p(ao = . — "m. — , (7) 

(s-M%) +M]p(s) 
where Mf denotes the fo mass and T (s) the width which is a function of the mass having 
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the explicit dependence 

r W =r„^EM.^, ,8, 

M% - AM 2 ) M ™ 

where M n is the charged pion mass, and T the width in the limit where the pion becomes 
massless. 

We find that M fo = (968 ± 9) MeV, and T = (92±|f) MeV. Later in this paper we will 
slightly refine these values of the mass and width, which we use for subsequent analysis. 
We note that performing this fit using a normal Breit-Wigner function gives Mf = (966 ± 
9) MeV, T = (89lljj) MeV. 



III. BRANCHING FRACTIONS AND q 2 DISTRIBUTIONS 
A. Main Experimental Results 

The main aims of this paper are to provide the branching fraction and q 2 dependence for 
the foe + v mode and to measure the relative rates at q 2 = of the /o and <fi modes. The 
BaBar collaboration previously measured B(D+ — > (j)e + u) = (2.61 ± 0.03 ± 0.08 ± 0.15)%, 2 
and the form factors for this decay given in the Appendix. They also found that the S-wave 
contribution to the K + K~e + u rate is a small fraction (0.221q;q|)% of the total rate in the 
mass interval between 1.01 and 1.03 GeV [3]. 

In general, for a decay of a pseudoscalar Df to a scalar meson such as the / (980), the 
only Lorentz invariants in the problem are the invariant square P 2 of the summed four- 
momentum of the D s and fo, and the invariant square q 2 of the four-momentum transfer 
between the D s and the f . The hadronic current describing the decay can be expressed as 

{fo | *y"(l - 7 5 )c | D s ) = U(Q 2 )P" + /-(?V , (9) 

where f±{q 2 ) are arbitrary form factors. Since this transition is between a pseudoscalar 
and a scalar f , only the axial- vector part of the current contributes. The term containing 
f-(q 2 ) gets multiplied by the positron mass squared, and becomes vanishingly small [14]. 
The decay rate can be written as 

dT (D+ -> f e + u) G 2 F \VJ 2 , . , 2 



dq 2 24tt 3 



Pfo 



U(q 2 ) , (io) 



where Gp is the Fermi constant, pf is the magnitude of the three-momentum of the /o in the 
rest frame, and \V CS \ is a Cabibbo-Kobayshi-Maskawa matrix element equal to ~ 0.97 
[1]. Therefore, we can use our data to determine |/ + (g 2 )|. 

Equation 10 is strictly correct only for zero width. It can be modified, however, to correct 
for the finite width mesons as described by Isgur et al. [15]. Here we integrate the product 
of P(s) as defined in Eq. 7 with the right-hand side of Eq. 10 over dM n7T . 



2 The first error is statistical, the second experimental systematic, and the third due to the absolute D+ 
branching fraction scale. 
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The form factors for 0e + z/ are somewhat more complicated as the hadronic transition 
involves a transformation from a pseudoscalar D s to a vector 0. The formulas for this case 
are given in the Appendix. The form factors for this decay have been measured rather 
precisely by BaBar [3], and we use these for further analysis of this mode. Our aim is to 
determine our own rate at q 2 equals zero, so as to cancel systematic errors in the ratio of 
rates for fQe + v/(f)e + v. 

We separate the data into five q 2 intervals and perform fits to the mass distributions. 
These are shown in Fig. 6 for 7r + 7r~e + z/, and in Fig. 7 for K + K~e + v. The resulting signal 
yields are given in Table III. For the 7i + n^e + u mode, the fitting function consists of a sum 
of a signal Breit-Wigner with mass and width fixed to the values we found by summing 
all q 2 intervals for the fo, a fixed background function from the D~ sidebands, a floating 
background shape to account for — > ?/e + z/, rj' — > 7r + 7r _ 7, and a fixed background function 
to account for — > 0e + z/, — > n + n^n . The K + K~e + v shape is defined in section II C 
and is kept the same in each q 2 interval. 



TABLE III: Number of events (#) and efficiencies e in q 2 intervals. 



q 2 interval #/o e+z/ 


<fo) 






<4>) 


#0e+^ 


(GeV 2 ) 




(%) 


corrected 




(%) corrected 


0-0.2 


14.6±3.9 


45.5 


32.1±8.6 


30.5+5.1 


25.9 


118±20 


0.2-0.4 


12.3±3.5 


50.2 


24.5±7.0 


30.5±5.3 


22.9 


133±23 


0.4-0.6 


12.4±3.5 


53.6 


23.1±6.5 


22.3±4.4 


21.5 


104±20 


0.6-0.8 


4.2±2.2 


56.9 


7.4±3.9 


21.5±4.7 


15.6 


138±30 


0.8-2.0 


0±1.0 


55.9 


0±1.8 


1.6±1.1 


9.8 


16±11 


Sum 


43.5±6.7 




87.1±13.5 106.4±9.8 




508±49 



The product branching fractions can be extracted by dividing the efficiency-corrected 
event sums in each mode by the number of D~ tags. We find 

B(D+ -> / (980)e + z/, f -> tt+tt") = (0.20 ± 0.03 ± 0.01)%, (11) 
B(Df^(f)e + u, <f)^K + K-) = (1.16 ± 0.11 ± 0.06)% , 

where the first error is statistical and the second systematic (see section III B). These values 
are consistent with, though the fo mode is somewhat larger than, our previously published 
values of (0.13±0.04±0.01)% 3 and (1.13±0.19±0.06)% for the /„ and modes, respectively. 
We can compare directly with the branching ratio measurement of the mode of BaBar [3] 
by dividing by the — > K + K~ branching fraction of (49.2±0.6)% [1]. Our measurement 
then becomes 

B(D+ -> 0e + z/) = (2.36 ± 0.23 ± 0.13)% , (12) 

while the BaBar measurement is (2.61 ± 0.03 ± 0.08 ± 0.15)%. 

Figure 8 shows the q 2 distributions for the foe + u and 0e + z/ channels. We fit the (pe + u 
channel using the form factors determined by the BaBar collaboration letting only the nor- 
malization float [3]. The fit yields 478±45 events compared with 508±49 from summing the 
data. 



3 Our result is larger here because the previous analysis assumed the fo natural width was 50 MeV. 
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FIG. 6: The invariant ir + ir~ mass in the semileptonic mode for the different q 2 intervals in units 
of GeV 2 : (a) 0-0.2, (b) 0.2-0.4, (c) 0.4-0.6, (d) 0.6-0.8 and (e) 0.8-2.0. The fits to the data are 
described in the text. The dotted curves show the signal, the long-dashed distributions are from 
sidebands of the D~ candidate invariant mass distributions, while the shorter dashed curves show 
the r( —* 7r + 7r~7 background level, and the dashed-dot curves show the (small) 7r + 7r~7r background 
from (be + v. The solid curves show the total. 



For fo£ + v we fit the q 2 distribution to Eq. 10 using a simple pole model for the shape of the 
form factor, |/ + (g 2 )| = 1/(1 — q 2 /M 2 olc ). (We exclude the highest q 2 point from the fit.) This 
procedure gives a value of (1.7±g;f) GeV for M pole . We then vary our measured fo width by 
one standard deviation and repeat the procedure to find the systematic uncertainty, resulting 
in a ±0.2 GeV systematic error. Using the central value of the pole mass we find that the 
fitted yield for foe + v is 86.6±13.0 corrected events compared with 87. 1±13.5 corrected events 
found by summing over the q 2 intervals. 

We now use these fits to determine Rf/<f,- At zero q 2 we find 38.8±9.3 corrected f§e + v 
events, and 92.6±8.7 <pe + v events. Thus 

R m = (42 ± 11)%. (13) 

Here any errors in the event tagging or reconstruction efficiencies cancel. The systematic 
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FIG. 7: The invariant K + K~ mass in the semileptonic mode for the different q 2 intervals in units 
of GeV 2 : (a) 0-0.2, (b) 0.2-0.4, (c) 0.4-0.6, (d) 0.6-0.8 and (e) 0.8-1.0. The solid curve show the 
total and the long-dashed distributions are from sidebands of the D~ candidate invariant mass 
distributions. 

error due to the uncertainties on the BaBar form factors are much smaller than the statistical 
errors. 



B. Systematic Errors 

The systematic errors on these branching fractions are given in Table IV. The error 
on track finding is determined from a detailed comparison of the simulation with double 
tag events where the entire event can be reconstructed even if one track is missed [16]. The 
particle identification uncertainty on the two hadronic tracks is listed as twice the uncertainty 
on the identification of one hadron. The error on the number of tags of ±2% is assigned 
by varying the fitting functions and ranges. Additional systematic errors arising from the 
background estimates are at the 1% level. Final state radiation effects are included in our 
simulations and the resulting accumulated uncertainty from the two-hadrons and the lepton 
is on the order of 1%. 
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> 




FIG. 8: The q 2 distributions for (a) Df — ► foe + u fit to a function allowing a varying pole mass 
and (b) Df — ► </>e + f fit to a function with form factors fixed to those measured by BaBar [3]. 



We have previously checked the resolution in MM 2 by fitting to a sample of Df — > K + K° 
candidates where we first found the K° and then ignored it [10]. The resolution agreed with 
the simulation to 1.7%. Here, because we have a positron in the final state instead of a 
K + we use an uncertainty of ±3.4% to allow for uncertainties associated with radiation in 
detector material. We have made wide selection cuts to minimize the effect of any resolution 
errors on our extracted rates. 

We note that there is no enhancement in our ability to find tags in /o e+z/ or 0e + z/ events 
(tag bias) as compared with events where the -D+ decays generically. Using a Monte Carlo 
simulation for each tag mode independently and then average the results based on the known 
tag fractions, we find a correction factor of (0±1%), which we assign as a systematic error. 



TABLE IV: Systematic errors on determination of the branching fractions. 



Error Source 


Size (%) 


Track finding 


2.1 


Hadron identification 


2.0 


Electron identification 


1.0 


MM 2 width 


3.4 


Background 


1.0 


Number of tags 


2.0 


Tag bias 


1.0 


Total 


5.2 
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IV. FURTHER RESULTS AND IMPLICATIONS 



A. Strange Mixing Angle of the / (980) 

In one class of models, the / (980) is thought to be a J p = + state described in the 
quark model as composed of quark-antiquark states which form a mixture of strange and non- 
strange components characterized by a mixing angle 9 defined in terms of the wave-function 
as described by Aliev and Savci [17] 

| f ) = cos 9 | ss) + sin#-^= | (uu + dcty) . (14) 

Since the central value of the mass is below threshold for decay into a kaon pair, the 
branching ratio to physical ss states is suppressed. BES has extracted relative branching 
ratios using ip(2S) — > 7x c o decays where the Xco —> fo (980) f (980), and either both / (980)'s 
decay into it + it~ or one into it + it~ and the other into K + K~ [18]. From their results we 
obtain 

Assuming that the tttt and KK decays are dominant we can also extract 

B(f - 7T+7T-) = (50^ 7 )% , (16) 

where we have assumed that the only other decays are to 7r°7r°, 1/2 of the 7r + 7r~ rate, and 
to neutral kaons, equal to charged kaons. This approach has been used by B. El-Bennich et 
al. [19]. It is possible, however, that the decays f — > 7r + 7r~7r + 7r~ or 7r + 7r~7r°7r° are significant. 
If so, then the branching fraction of fo into two pions will be smaller, 

Using the above branching fraction we extract B(Df — > foe + u) = (0.40 ± 0.06 ± 0.06)%, 
where the error on the fo — > 7r + 7r _ rate has been included in the systematic error, and, in 
fact, is the dominant contribution. 

We can use this rate to estimate the mixing angle cos 9. Aliev and Savci using QCD sum 
rules predict that B(D+ — > foe + u) = cos 2 9 x (0.41)%. This provides us with a value of 

cos 2 9 = 0.981°;^ , (17) 

where we have truncated the positive uncertainly so that the value does not exceed unity. 
This determination gives a somewhat larger but consistent evaluation with other methods 
[19]. If the fo decay rate into four pions is significant, our value of the mixing angle will be 
somewhat larger. 

The near unity value for cos 2 9 poses a conundrum: If the fo is dominantly an | ss) state, 
why does it decay predominantly into ttttI A solution may be that the quark content of 
the fo is not composed of just a quark-antiquark state, but may contain a some four-quark 
content [2]. 



B. The / (980) Mass and Width 

We now return to the determination of the fo mass and width. In the fit described above 
we used a relativistic Breit-Wigner function. We did not, however, take into account that 
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the phase space for the larger f masses is somewhat smaller, due to the finite D s mass, 
than for smaller / masses. The changes to the invariant mass spectrum of the /o have been 
parameterized by Dosch et al. [13]. The double differential decay rate in terms of s = M% n 
is 

T^iiS: A3/2(M -- s -' 2)l/+( ' 2)|2pw - (18) 

where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz and is equal to 4M|, a times the square of the 
momentum of the fo in the D s rest frame, and P(s) is the relativistic Breit-Wigner shape 
given in equations 7 and 8. The resulting, slightly shifted mass spectrum, then is given by 
integrating the above expression over q 2 

_*L = 2V~s / (MDs "^ )2 f^l d q 2 (19) 
dyfs Jo dsdq 2 

We perform this integral using our fitted form factor shape for |/ + (g 2 )| and our original 
fitted parameters for the relativistic Breit-Wigner function. This generates a new shape 
for the mass spectrum, that we fit to the data. In principle this procedure can be iterated 
many times; here we find that only one iteration is sufficient. The resulting curve fitted to 
the data is shown in Fig. 9. The systematic errors on the mass are small compared to the 
statistical error. The masses of known particles, K$, A, 0, D° are reproduced in other CLEO 
analyses to better than 1 MeV. Changing the background level by ±1 standard deviation 
changes the mass by a negligible amount. For the width we have identified two sources of 
systematic error, one the resonance shape and the other the form factor shape. Changing 
the Breit-Wigner shape from relativistic to non-relativistic causes a change of ±2 MeV. We 
also vary the form factor by ± one standard deviation resulting in a change of about ±2 
MeV. Our final values for the mass and width are 

M fo = (977+1 1 ± 1) MeV, and T = (91^ ± 3) MeV . (20) 

We do not fit our data to a function that allows for the opening of a threshold in another 
channel, because there is no obvious change in the shape of the dipion mass spectrum near 
1 GeV, and we do not have enough data to specify the parameters of such a fit. 



V. CONCLUSIONS 

We present an updated result for the first measurement of 

B(Df -> f (980)e + u)B(f -> n + n') = (0.20 ± 0.03 ± 0.01)%. (21) 

Assuming a simple pole model for the form factor, |/ + (g 2 )|, we estimate the pole mass as 
(1.7^;f ± 0.2) GeV. 

In the final state the only hadron present is the /q. This provides a particulary clean 
environment that allows us to measure the mass and width as (977lg 1 ±l) MeV, and (91 ji^ =t 
3) MeV, respectively. 

We update our measurement of 

B(D+ -> (j)e + u) = (2.36 ± 0.23 ± 0.13)% . (22) 
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n + n~ Mass (GeV) 

FIG. 9: Invariant mass distribution for 7r + 7r~ in the semileptonic mode. The dotted curve shows 
the signal shape, a relativistic Breit-Wigner modified by phase space and form factor effects, the 
long-dashed distributions are from sidebands of the D~ candidate invariant mass distributions while 
the shorter dashed curve in shows the rf — * 7r + 7r~7 background level, and the dashed-dot curve 
shows the (small) 7r + 7r~7r° background from (pe + v. The solid curve shows the total. 



We measure the ratio of decay rates at q 2 = 



dr 

„2 



(£>+ - /o(980)e+z/, f - tt+tT 



R fl* = 7^775 . 77777 1 . r/ + r/-V ( M ) 



This ratio has been predicted by Stone and Zhang to equal [6] 

T(£? s -> J/^/ , fp -> 7T + 7T- 

Our measurement indicates that the 5 S — > J/ipfo channel may indeed be a useful place to 
measure CP violation in the B s system in that the rate can be ~40% that for J/i/xj) mode, 
especially since an angular analysis is not necessary, as the J/V'/o mode is a CP eigenstate. 

Finally, we show that viewing the fo as a quark-antiquark state results in a mixing angle 
between ss and light quarks whose cosine is close to unity, consistent with other estimates 
[19]. Thus, the large rate into tctt may be the result of part of the wave-function being a 
four-quark state [2]. 
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APPENDIX 

We briefly summarize the relationships between different rates and form factors in D s 
decays into a vector particle (0) here. There are more Lorentz invariants possible than 
in the pseudoscalar to scalar case, and three independent form factors, V, A±, and A2 are 
involved in the case of e + u decays. The matrix element can be written as the sum of a vector 
plus axial- vector current [14]: 

V(a 2 ) 

(<t>\SlADs) = 2 Md Im/ ^PDsPV^ and (25) 
(01^75^) = i (M Ds + M+) {el - ^-q^j A^q 2 ) (26) 



M Ds + M 4> \ q 2 

where e* is a polarization of the final state meson. (An additional term proportional to Aq 
has been dropped as it gets multiplied by the square of the electron mass.) 

First of all the differential decay rate can be written in terms of the helicity amplitudes 
as [20] 



dT G 2 F \V CS \ 2 



dq 2 96tt 3 M 2 Dt 



(|# + | 2 + |#_| 2 +|# | 2 ) • (27) 



The helicity amplitudes are expressed in terms of three form factors in the limit of zero 
lepton mass as [21] 



H (q 2 ) ' 



2M D 



(M 2 Ds - M 2 - q 2 ) (M Ds + M^q 2 ) - * Md °+ m M<?) 



2 ^2 

' , (28) 



H ± (q 2 ) = (M Ds + M^Mq 2 ) T (Md +m ) y (g ) ■ ( 29 ) 

The most elementary parameterization of the form factors considers the q 2 dependence 
as the behavior of a simple pole in the t-channel, yielding 

The masses can be taken as My = Mr,*, and Ma = Md s1 , although they can, in principle, 
be extracted from the data. The BaBar collaboration has assumed the value for My and 
measured M A = (2.281q; 2 | ± 0.18) GeV; they also determined at q 2 = 0, r v = V(0)/Ai(0) = 
1.849 ± 0.060 ± 0.095; r 2 = A 2 (0)/A 1 (0) = 0.763 ± 0.071 ± 0.065 [3]. 
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